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^ ' Abstract 

We discuss the BFKL approach to processes with large momentum transferred through 
a rapidity gap. The Mueller and Tang scheme to the BFKL non- forward parton-parton 
elastic scattering amplitude at large t, is extended to include higher conformal spins. The 
new contributions are found to decrease with increasing energy, as follows from the gluon 
reggeisation phenomenon, and to vanish for asymptotically high energies. However, at 
moderate energies and high the higher conformal spins dominate the amplitude. We 
illustrate the effects by studying the production of two high Ej- jets separated by a rapid- 
ity gap at HERA energies. In a simplified framework, we find excellent agreement with 
the HERA photoproduction data once we incorporate the rapidity gap survival probabil- 
ity against soft rescattering effects. We emphasize that measurements of the analogous 
process in electroproduction may probe different summations over conformal spins. 



1 Introduction 



High energy hadronic scattering with large momentum transfer |t| ^ ^qcd ^^"^ rapidity y ^ 1 
is an excellent testing ground for perturbative QCD. The most interesting is the case when 
colour is not exchanged in the interaction. Our understanding of such processes (hard colour 
singlet exchange) is based on the BFKL equation |I|, 0] which resums the gluonic ladder 
diagrams in the leading logarithmic approximation^. The best known processes in which the 
behaviour of the scattering amplitude may be tested are elastic vector meson production 
^ 0, diffractive, proton dissociative •yp scattering p and events with gaps between jets in 
the appropriate kinematic regime |T^, [T^, |T^. Two main approaches to determine non- 
forward BFKL amplitudes have been proposed. 

The first approach relies on the conformal symmetry of the leading logarithmic BFKL 
equation, which permits an analytical solution of the problem WM, R[. The applications follow 



the Mueller and Tang |^ subtraction scheme to obtain the elastic parton-parton scattering 
amplitude. It is valid for an asymptotically large rapidity gap. However, it is doubtful whether 
this asymptotic formula may be used for the currently available measurements. The main 
problem is, that the Mueller- Tang cross-section for gaps between jets is, even for ?/ ~ 5, much 
smaller than the lowest order two-gluon exchange cross-section. Of course, one has to include 
the gluon reggeization factor which suppresses the infrared sensitive part of the two-gluon 



amplitude [H, Q5|, but still this contribution appears to be very important |12, 13] for the region 



of y probed in the current experiments [ p^ . 

The other approach is based on numerical studies of the non-forward BFKL equation 0, 



1^]. The main advantage of this method is that it does not require any restrictions imposed on y 
and takes into account all the available details of the impact factors. An important ingredient 
of this framework is that it is possible to go beyond the leading logarithmic approximation 
by including some phenomenological modifications of the BFKL kernel which are expected to 
resum a major part of the higher order corrections, like the running of the coupling constant 
along the ladder and the imposition of the consistency constraint [|1^, 0, 0] . 

The main purpose of this paper is to generalize the analytical approach so that the BFKL 
parton scattering elastic amphtude can be used at lower values of y. We will demonstrate 
that in this case we have to include the higher conformal spin contributions and so we are 
able to trace the phenomenon of gluon reggeization in a representation given by the conformal 



^Recently, also the next-to- leading corrections have become available 
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Figure 1: The Feynman diagram illustrating the mechanism of photo- and electroproduction of two 
high Et jets separated by a rapidity gap. 

eigenfunctions of the BFKL kernel. Possible phenomenological effects of higher conformal spins 
in the forward BFKL amplitude have been discussed, for example, in [|^, 

From the experimental point of view, the process of interest (with a large rapidity gap 
between two high Et jets) has been observed at the Tevatron and now data are becoming 



available for the analogous process at HERA ETI, |22[. For example, in Fig. 1 we show the 



production mechanism for the diffractive process at HERA. Originally, it was claimed [0 that 
the Tevatron observations strongly disagree with the BFKL approach. However a closer study 
shows that the contradiction disappears when we allow for the effects of (i) hadronisation, 
(ii) the survival probability of the rapidity gaps and (iii) asymmetric non-asymptotic BFKL 



contributions. Numerically, at these energies it was found at the partonic level [p!2| , |T3[ that the 
elastic parton-parton amplitude may be well described by two reggeized gluon exchange. Indeed, 
as we will show in Section 2, the non-asymptotic (i.e. higher conformal spin) components may 
be summed to give reggeized two gluon exchange. We note that the conformal spin components 



contribute up to the value of the conformal spin n ~ ^ET/ko, where the infrared physical cutoff 
ko is driven by the size of the incoming state. 

In Section 3 we consider the corresponding process at HERA: that is two high Et jets 
separated by a rapidity gap as shown in Fig. 1. Here we have a second variable, the photon 
virtuality Q^, which allows us to change the size of the incoming qq state. In this way one can 
study the effect of varying the number of conformal spin components. We have attempted to 
summarize the essential points of the analysis in a self-contained, and more physical, way when 
describing the application in Section 3. 



2 Elastic parton-parton scattering amplitude 



The quark-quark elastic cross-section at high energies reads 



'^'''Hy) = ^\A{y,t)\^ (1) 



dt ^ ' SItt 
where the amphtude Aiy^t) is 

A{y,t) = I d'kd'k'r{k,k';y). (2) 
2.1 Two gluon exchange 

In the ?/ — > hmit, the BFKL amphtude reduces to its lowest level approximation of the 
exchange of two elementary gluons, i.e. 

The integral over the tranverse momenta in eq. (|^) is infrared divergent in this limit. One may, 
however, introduce an infrared cutoff k^ and modify the gluon propagators + ^o)- 

Such cutoff has a physical origin. It may be related either to hadron sizes or to the gluon 
propagation length in the QCD vacuum. With such a substitution, becomes 

27r 

A{y = 0,t)^-\og{qykl) (4) 

where q is the momentum transfer. For the production of a pair of high Et jets, we have 
q ~ Et- At not large rapidities and large \t\, the most important effect of the BFKL evolution 
is that the exchanged gluons become reggeized. This accounts for a no-emission amplitude 
from a system of two gluons in the colour singlet state, where one of the gluons is much softer 
than the other. To a good approximation, it leads to an additional factor multiplying the gluon 
propagator close to the singular point 

1 1 



(5) 



with z = SasV/ (Stt). Then for moderate y 

X 27r fi' „2 1 fk^Y 

For y ^ the integration in (M) may be safely performed, and the limit k^^O taken, to obtain 



X 27rl 
A{y,t = — - 
q^ z 



% (7) 

zq^ 



We shall trace how the above expression arises from the summation over the conformal spins. 
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2.2 BFKL conformal components 



The LO BFKL expression for p{k, k'; y) is |T|, § 



X 



(2vr)6 J {[i^' + in- 1)74] [z.^ + (n + 1)74] 

exp[^„(z/)y]/i,,(fc,g)/2/(fc',q)} (8) 

where 

Un{jy) = —[2^(1) - ^(1/2 + \n\/2 + iv) - ^(1/2 + \n\l2 - iz/)] (9) 
7r 

are the eigenvalues of the BFKL kernel. The functions ^ are constructed from the impact 
factors q) and the eigenfunctions of the BFKL kernel. To be precise 

/^^^(fc, q) = q) J (fpi (fp2 EnApuP2) exp{ikpi + i{q - k)p2) (10) 

where the eigenfunctions take the form 

E„AP.P.)=(^)''(i^)'t (11) 

V PlP2 J \\ Plp2 J J 

Here we have used the complex representation of transverse vectors k and p (which reveals the 
conformal symmetry), that is k = + iky and p = Px + ipy The powers h = 1/2 + n/2 + iu 
and h = l/2 — n/2 + iv are the conformal weights. 

For parton scattering, the impact factors $*(A;, g) are constant functions of k. In other 
words for a point-like parton we have pi = p2, and hence the eigenfunctions (|TTD , the impact 
factors (|l^) and the amplitude of (^, vanish identically. However one can not use the ex- 
pansion (^) over the conformal eigenfunctions En^u for coloured initial objects. In such a case 
we face an infrared divergency in the integrations over the pi which correspond to contributions 
proportional to 6{k) or 6{q — k). For a colourless object these divergent terms are absent, 
since they are multiplied by zero - the total colour charge. Thus one has to consider the par- 
ton spectators which compensate the colour charge of our active quark. Assuming that these 
spectators are located at rather large distances we may follow the Mueller- Tang prescription. 



as was shown in 15 



The generalisation of Mueller- Tang prescription for arbitrary conformal spins reads: 

'pi-P2y((pi-P2\y fpi-p^yffpi-pAy fMVM' f-^Vf-^V 

. P1P2 J \\ Plp2 J J ^\ P1P2 J \\ Plp2 J J \P2J [py \PlJ \PlJ 
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, Plp2 I \\ Plp2 I I \\P2\I \\P2\I \\pl\) VIPiI/ 



(12) 



Note the minus sign in the last term, which will result in the cancellation of contributions with 
odd n. We set fc) = 1 for s = 1, 2 and substitute (|12D into ([T0|) to obtain 



We insert the last result into eq. and then integrate over k and k' to obtain the non-forward 
amplitude: 



4 °° r ( 

Aiy^i) = -2 E /^H 

y m=-oo'' I 



[z/2 + (m - 1 /2)2] [zy2 + (m + 1 /2)2] 



exp[u;2m(«^)2/] [ • (14) 



The last formula represents the desired generalisation of the Mueller- Tang result for the quark- 
quark elastic scattering amplitude. Note, that only contributions from even conformal spins 
n = 2m are left in the sum. 



2.3 Sum over conformal spins 



The resulting amplitude has the following properties. For very large y all the components with 
|m| > get suppressed because in this case uj2mW) < 0. Then, indeed, it is enough to retain 
only the leading term with m = (i.e. n = 0), which gives rise to an increasing part of the 
amplitude with the famous LO BFKL intercept 



A{y,t) ~ y ^'/^exp 



12 log 2 as 



TT 



■y 



(15) 



However, for ?/ — > the expression is divergent, as expected from an inspection of the two gluon 
exchange amplitude, due to high l/\m\ asymptotics in the terms under the sum (H): 



dv 



z/2 + m? 



n 



TT — 1 / J 



for 171 = 0, 
for m 7^ 0. 



(16) 



+ (m - 1 /2)2] [z/2 + (m + 1 /2)2] 

|m| m2 — 1/4 

For y > the divergence disappears due to the presence of the exp[co'2m('^) y] term. Namely, for 
large m one has 



UJ2m{l^) < ^2m.{l' = 0) 



6a.o 



TT 



log(|m|)-7E] + 0(l/m2), 



(17) 
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Figure 2: The integrals over v in ^1^ ) corresponding to values of m = \n\/2 between and 10 (dashed 
lines), the sum over the integrals from n = —4 to 4 and from n = —20 to 20 (continous lines) and the 
naive sum of n = component and the reggeization term (Qj (its l.h.s.) multiplied by g^/4, plotted as 
functions ofSasU/T^ in three cases: q/ko = 4? , 20^ and ko = (dotted lines). 



which bounds from above the suppression factor in the sum over m to be |m|-6^°=/^. This 
changes the l/\m\ asymptotics of the summed terms and ensures the convergence of the infinite 
sum. Note, that for small y (i.e. z —* 0) the limiting behaviour of the amplitude is 

oo StT 

A{y, t) ~ — m^"'^"''^ = — - + regular terms. (18) 
m=i T zq' 

Thus we have reproduced the answer (^,^, which was obtained by accounting for the gluon 
reggeization when the infrared cutoff ko 0. 

In Fig. 2 we illustrate the result of taking a finite number of terms in the sum over the 
conformal moments in ([T^) . We first plot the sum of terms with all conformal spins n between 
-4 and 4, and then show the sum for terms between -20 and 20. Only even values of conformal 
spins contribute. 
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2.4 Insight from comparison with naive estimates 



At very high rapidity the amphtude is dominated by the n = component, while at low rapidity 
it is mainly two gluon exchange (^. We may therefore try to approximate the full amplitude by 
the sum of these two contributions. To compare this approximation with amplitude summed 
over conformal spins, we have to understand how the infrared parameter reveals itself in the 
BFKL expression which was written in terms of conformal eigenf unctions. It turns out that 
the impact factor In,u goes to zero for \n\ > rimax, where the value of nmax is regulated by the 
ratio q/kQ. 



Let us focus on the integral ([T0|) defining the function ^(fc, q) for configurations with, say, 
k <^ q, which dominate at low y. In this case the typical values of pi and p2 in integral (p!0|) 
are |pi| ^ |p2|. The scale for |pi| is set by the size R of the initial colour multipole. 

Then, one may rewrite 

V pip2 J VV piP2 / / \p2j \pi) V pi pV 

where only the leading correction in |pi/p2| to the result obtained in the limit R ^ oo, is 
retained. Thus, a rough estimate of the relative correction to (|1^) coming from the fact that 
the size R is finite, from the term proportional to hp2/pi , is given by the ratio 

^ _ -Id^p2 {hp2/R) exp(igp2) Pa^Pa)"^ ^20) 

/rf2p2 exp(igp2) P2''(P2)"'' 

The integrals in ( pUD are of the same form as those used for the derivation of (|T3[). In particular, 
the result for the denominator is already known to be 

27ri" /g\2i^r(l/2 + n/2-iz/) 
q Vlj T{l/2 + n/2 + iu) ^ ' 

and the integral in the numerator may be obtained from (^Tj) by substitutions iv iv — 1/2 
and n ^ n — 1. Thus, it is straightforward to find that 

iqR 

It may be seen that the typical value of \v\ in integral (H) is \v\ ~ \n\/2 for large \n\. Then 
\{n/2 + iz/)2 — 1/4| ^ 71^/2 and after including a similar result from the term in eq. ([T9|) 
containing hp2/pl , the estimate of the total correction is 

2\CnA - -5- (23) 
qR 



The Mueller- Tang scheme breaks down when the correction factor becomes large 2\Cn,v\ ^ 1 
and J^j, gets suppressed in relation to the r.h.s of (|T3|). Noting that ~ 1/-R^, we find from 
I) that the sum over conformal spins should be extended to 

/ q 

^max = \ 1~- (24) 

V 

We may compare the result of the summation up to rimax with naive estimates given by 
a simple addition of the n = contribution and the amplitude given by the two reggeized 
gluon exchange (|). The dashed curves in Fig. 2 show results for three values of the ratio g/fco, 
namely 4^, 20^ and oo. The agreement between the naive estimates and the truncated sum 
over conformal spins can be understood as follows. At y = the sum over n in (0) takes form 



A(y = 0,t)^^^U+ V 3 -^log , (25) 




where we have used (JT^ and (]TB|). In this way we reproduce the leading behaviour of the 



amplitude given by @. Indeed, as seen in Fig. 2, the naive sum of the n = component and 
the two (reggeized) gluon contribution (|^) for the two cases of q/k^ = 4^ and 20^ are rather 
close to the sum over the conformal spins \n\ up to 4 and 20 respectively^. 

Note that we do not have any non-perturbative effects modifying the gluon perturbative 
propagators. However, the sensitivity to the infrared details decreases with increasing y, since 
these details influence more the contributions from higher conformal spins, which get suppressed 
with y. The suppression is faster for larger \n\. In view of the above discussion and the recent 



results of Ref. ||T2|, |T3|, the extension to include higher conformal spins is necessary to understand 



the Tevatron data for hard colour singlet exchange. 



3 Phenomenological consequences for gaps between jets 

As we have just discussed, the above formalism is relevant to processes mediated by colour 
singlet exchange with high momentum transfer q. As noted, the classic example is the pro- 
duction of a pair of high Et jets separated by large rapidity gap y. We emphasize that the 
conventional (asymptotic) BFKL amplitude only dominates at high rapidities well above the 

^Note that, in general, the naive sum sHghtly overestimates the true result as the two gluon contribution 
already contains part of the n — component. 
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reach of present experiments. This is illustrated by the discrepancy between the n = curve 
and the full result (given by the continous curves on Fig. 2), in the region "iasy/ix < 1 currently 
sampled experimentally. Instead, in this region it is necessary to sum all the contributions with 



conformal spins up to n^^^ = ^ q/k^, where fco is a physical infrared cutoff. On the other hand, 
in this domain the amplitude is well described by two reggeized gluon exchange, as illustrated 
by the dotted curves in Fig. 2 (note, that the uppermost curve corresponds q/ko ^ oo whereas 



the lower dotted curves correspond to the choices yq/k^ = 4 and 20). The physical meaning 
of gluon reggeization, that is of the factor {kP' /q^Y in (|^,H), is that it reflects the fact that the 
emissions of extra gluons are forbidden within the rapidity gap interval so that pure elastic 
parton-parton scattering occurs. This is equivalent to the normal Sudakov-like suppression. 
The latter suppression is the probability not to emit gluons with transverse momentum pt in 
the interval {ko,q), which takes the form exp(— n^). The quantity n^, the anticipated average 
number of emissions, is 

^P? 3a, 

^9= L ~ — y- (26) 

Thus, the non-forward lowest-order two-gluon exchange amplitude should be multiplied by the 
suppression factor 

exp(-n,/2) = (^-f j , (27) 

as given in (|^). Thus we see that in the BFKL amplitude the suppression is generated by the 
resummation of the virtual corrections. 



The hard colour singlet exchange has been investigated experimentally at the Tevatron Wt 



by measuring events with gaps between jets. In [O] these data were compared with BFKL re- 
sults in the standard Mueller- Tang approximation. It was found that the rising dependence 
of the gap fractions may be reproduced by the model only when a fixed value of coupling con- 
stant is used. The assumed lack of running of the coupling with increasing momentum transfer 
is, however, difficult to motivate. The experimental status of the gap fraction dependence on 
the jet separation in rapidity. Ay, is not so certain as the observed Et dependence. In partic- 
ular, CDF results indicate a decreasing tendency at high Ay, contrary to DO data. The error 
bars are still too large to claim inconsistency and both the experimental distributions agree 
with a flat Ay dependence. This is, however, incompatible with the predictions based on the 
Mueller- Tang approximation, which give a steep rise of the gap fraction at large rapidity |TT 



In the Tevatron kinematical conditions the bulk of data comes from Et ~ 20 GeV and 
Ay ~ 5. Then, if our choice of the cutoff scale /cq ~ 1 GeV^ is correct, the Mueller- Tang 
approximation gives about a half of the contribution to the scattering amplitude, thus about 
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Figure 3: Gap fraction for photoproduction of two high Et jets separated by a rapidity gap Arj. The 
experimental data are obtained with the rapidity gap events defined as those with the maximal total 
transverse energy flow in between the high-Ex jets to be E^^ = 0.5 GeV. The continous and dashed 
curves show the partonic non-leading BFKL prediction, with and without the gap survival probability 
factor included respectively, and the dotted curve is the LO BFKL prediction without the gap survival 
probability factor taken into account. 



25% of the cross- sect ion. After setting as = 0.17, as in [|rT|, this may be seen in Fig. 2 by 



comparing the sum over conformal spins up to n = 4 with n = component, at 3asy/7f ~ 0.8. 
Therefore the Mueller- Tang approximation should not be used to describe the available Tevatron 
data unless the cut-off ko is much larger than we expect. 

Recently |T2|, [T^ it has been demonstrated that when the gluon reggeization phenomenon is 



accounted for, no significant discrepancies between the data and BFKL results appear in either 
the Et or the Ay distribution. This conclusion holds both for fixed and running couplings, 
provided they are consistently used in the amplitude at the scale set by the typical virtuality 
of each vertex. 

The result obtained in Refs. [O, 0^, from theoretical considerations based on the BFKL 



framework, was confronted with experimental data from the CDF and DO collaborations |T6 
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Besides the hard parton-parton scattering amphtude for colour singlet exchange, other im- 
portant effects, like hadronisation corrections and the gap survival probability were included 
using a complete Monte Carlo treatment. It was found |T^, [T^ that the full BFKL prediction 



was in agreement with the Tevatron data, contrary to calculations based on the asymptotic 
Mueller- Tang approximation. 

An interesting way of studying this effect in more detail is to observe events with a large 
rapidity gap between two high Ex jets in diffractive deep inelastic scattering at HERA, as 
sketched in Fig. 1. Since one jet comes from photon dissociation, the infrared cutoff ko is 
controlled by the photon virtuality Q^. Therefore, there is the possibility to vary ko while 
retaining the same kinematics of the hard parton-parton interaction. In this way, we avoid 
complications from hadronisation, variation of parton densities etc. Another advantage of 
electroproduction is that we have high survival probability, S*^, of the rapidity gap against soft 
rescatterings, that is S*^ ~ 1, contrary to the analogous jet production process at the Tevatron. 

At present, however, data are available at HERA for the photoproduction of jets sepa- 



rated by a rapidity gap. These data are compared with the corresponding non-forward BFKL 
predictions in Fig. 3. The two uppermost curves are obtained using the solutions to the BFKL 
equation given in Ref. |T^, and correspond to BFKL amplitudes with and without incor- 
porating resummations of higher order effects. In the leading order BFKL calculation a fixed 
value of as = 0.17 was used, whereas the running coupling was taken in the non-leading BFKL 
case, which explains why the LO BFKL curve lies so low. Nevertheless, for photoproduction the 
probability of soft rescattering, which produces secondaries in the rapidity gap, is not negligi- 



ble. To calculate the resulting suppression factor S we have used the formalism of Ref. ||23 
Recall that, there, the model was tuned to describe the available soft pp and pp interactions 
throughout the CERN ISR-Tevatron energy range. Assuming Vector Meson Dominance, the 
rescattering in photoproduction may occur between a virtual vector meson and the proton. 
Data are not available for soft vector meson-proton scattering. However using the additive 
quark model, and the analogy between the light vector mesons V = p,u!, ... and the pion, we 
expect^ 

crtot(^p) ^30mb (28) 
at the collision energies relevant to the HERA data, that is ~ 200 GeV. With this cross 



•^Note that the ZEUS coUaboration |2^] estimate that (Ttot(7i'p) = 31 ± 4 mb, by observing the interference 
between the pions from p meson and tt+tt^ background. 
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section we determine the relevant photoproduction rapidity gap survival probability to be [23 f\ 



~ 0.3. (29) 

After taking this factor into account we obtain the final prediction shown by the lower contin- 
uous curve in Fig. 3. There is thus an excellent agreement between the data and the theory. 

It is appropriate to comment on the approximations made to obtain this prediction for the 
photoproduction of jets separated by a rapidity gap at HERA. We have assumed that the effects 
of hadronization and of producing gaps in the conventional colour-octet exchange scattering 
are small. We expect the experimental cut on soft secondaries to suppress these effects. 

Moreover, although we believe our two-channel eikonal calculation of the survival factor 5*^, 
using the framework of [^|, is the best that can be done at present, we note that again it relies 



on soft phenomena. Clearly it is important to study all these effects in more detail within a 
Monte Carlo framework as in [0, |13|. This will, among other things, allow a study of the 
influence of the parameter i?™* in the gap definition on the gap fraction. 

When the HERA luminosity increases, it will be particularly informative to observe elec- 
troproduction of high Et jets separated by large rapidity gaps. This will allow Q^, as well 
as the jet Et-, to be varied, and hence conformal spin summations up to different n^ax to be 
probed. Moreover, here, there are no uncertainties connected with the survival factor S*^, since 
it is predicted to be 1 for this process. These data will therefore allow a test of QCD radiative 
effects, which are important ingredients in all non-forward diffractive phenomena. 
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